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Abstract 

When the electron density is a linear function of distance, it is known that the 
MSW equations for two neutrino species can be solved in terms of known functions. 
It is shown here that more generally, for any number of neutrino species, these MSW 
equations can be solved exactly in terms of single integrals. While these integrals 
cannot be expressed in terms of known functions, some of their simple properties are 
obtained. Application to the solar neutrino problem is briefly discussed. 
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1 Introduction 



In studying the Mikheyev-Smirnov-Wolfenstein (MSW) effect |J due to the coherent for- 
ward scattering of neutrinos by electrons in matter, it is often instructive to consider first 
special cases where the electron density is taken to be a simple function of distance. It is 
the purpose of the present paper to investigate perhaps the simplest case: the case where 
the electron density is a linear function of distance. 

The problem of the linear electron density is formulated in Sec. 2. The case of two 
neutrino species has a long history ||, and the solution, as reviewed in Sec. 3, can be 
expressed in terms of parabolic cylinder functions — see, for example, Chapter VIII of ||, 
or equivalently confluent hypergeometric functions. However, the solution in this form is 
specific to the case of two neutrino species, and is not convenient for generalizations to 
more neutrino species. Physically, this generalization is essential because there are at least 
three types of neutrinos. Therefore, Sec. 4 is devoted to treating in a different way the 
MSW differential equations for linear electron density and two neutrino species. On the 
one hand, this alternative method must lead to the same solutions as those in Sec. 3; on 
the other hand, this new treatment can be readily generalized to any number of neutrino 
species. This case of linear electron density but any number of neutrino species forms the 
main content of the present paper, and various aspects of this case are treated in Sec. 5 
and Sec. 6. 

2 Formulation of the Problem 

Let there be iV types of neutrinos, denoted by U\, . . . v N , where v x is the neutrino 
of the first generation, i.e., the one that forms the SU(2) doublet with the electron. It is 
assumed that v\ is the only neutrino which interacts differently with the electron because 
of the exchange of the intermediate boson W, while the others neutrinos i>2, ^3, • • • vn all 
have the same interaction with the electron. Thus, the neutrino mass matrix M 0] is an 
N x N matrix. The eigenvalues of M give the N neutrino masses fi. 

In analyzing the MSW effect, the neutrino masses are usually taken to be much smaller 
than the momentum p of the neutrino. Under this assumption, because 



it is M 2 that enters in the differential equation for the MSW effect. Let ^>(x) be the 
^-component neutrino wave function, then this differential equation is 



(p 2 + ^)V 2 ~p + -L^, 



(2.1) 



i— = W(x) + ^-M 2 



dx I 2p 



(2.2) 



where W(x) is an N x N matrix whose only non-zero element is 

[W(x)] n = V2G F N e (x), 



(2.3) 
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with Gf the Fermi weak- interaction constant and N e (x) the electron density at the point 



x. 



The terminology "linear electron density" is used to mean that N e (x) is a linear function 
of x. 

Since N e (x) is the density of electrons, it cannot be negative. Therefore, the MSW 
differential equation ( |2.2| ) is physically meaningful only for the half-line of x where N e (x) > 
0. On the other hand, when the neutrino or the electron, but not both, is replaced by its 
antiparticle, the quantity [W(a;)] 1:l of Eq. (|2.3j) changes sign. Therefore, the complementary 
half-line of x describes this slightly different physical situation. For this reason, Eq. (|2.2j ) 
is to be studied for the entire range of x, from — oo to +00. 

For the present case of the linear electron density, Eq. (|2.2j ) can be reduced, for a given 
value of p, to the dimensionless standard form 



ij t m = A{t)i>{x) 



where 



(2.4) 



(2.5) 



and 



A(t) 
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(2.6) 



This is accomplished as follows. 

(i) To change the independent variable from x to t, there is a shift in origin and a 
rescaling with possibly a reversal of sign. 

(ii) To change the dependent variable from \I/ to ip, there is a rotation in the second 
to iVth component and an introduction of exponential factors with possibly some minus 
signs. 

Furthermore, from (i) and (ii), the elements of the matrix A(t) can be chosen to satisfy 
the conditions 



TV 

i=2 



(2.7) 
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b 2 < h < W < ■ ■ ■ < b N -! < b N , (2.8) 

and 

dj>0 for j = 2,3,4,... N. (2.9) 
Consider the following special cases. 

(a) If, for some j, say j , a J0 = 0, then it is seen from Eqs. ( |2.5| ) and ( |2.6| ) that ipj * s 
decoupled from the other ^ 's. Thus, this special case of iV types of neutrinos is reduced 
to a problem with N — 1 types of neutrinos. 

(b) If, again for some j, say jo, bj = bj 0+ x, then a rotation can be carried out between 
ipj Q and ipj +i such that, after this additional rotation, the new a JO is zero. Thus, this 
special case of bj = bj 0+ i can be reduced to the above case of a 3 - Q = 0, and hence again 
this second special case of iV types of neutrinos is reduced to a problem with N — 1 types 
of neutrinos. 

It is therefore sufficient to study the ordinary differential equation (|2.4j ) with Eqs. (|2.5| ) 
and (|2.6|) under the condition (|2~7|) together with 



b 2 < b 3 < 6 4 < . . . < &7v_i < b N (2.10) 

and 

aj>0 for j = 2,3,4,... N. (2.11) 
In view of the inequality fl2.10p , it turns out to be convenient to define symbolically 

bi = —00 and b N+ i = +00. (2-12) 



3 Case N = 2 



Let us review first the well-known case of the MSW effect for two types of neutrinos [|T], . 
By Eqs. (O)-(O), the MSW equations are 



(3.1) 



d 






-t a 2 






I— 

dt 






a 2 







or more explicitly 



. d , , . 

, Tt Mt) 



-ty x {t) + a 2 ip 2 (t), 



(3.2) 
(3.3) 



with a 2 > 0. A second-order ordinary differential equation for ipi(t) is obtained by applying 
d/dt to Eq. Q and using Eq. (O) : 



d^tit) #iW 2 
- it — V (a 2 



dt 2 



dt 



i)Vi(*) = 0. 



(3.4) 
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In order to remove the first- derivative term, let 

ij 1 (t) = e u2 / 4 Mt) (3.5) 

Then the equation for <j)i(t) is 

^^ + (^ 2 + «i-^)0iW = O (3.6) 



Two linearly independent solutions of this Eq. ( |3.6j ) are the parabolic cylinder functions 

D p {±e**' 4 t), (3.7) 

where 

p = -ia 2 2 - 1. (3.8) 

Parabolic cylinder functions are special cases of the confluent hypergeometric function 
0, the relation being 

D p (z) = 2^)/ 2 e^ 4 z *(I - \p, |; ±z 2 ). (3.9) 

Since the confluent hypergeometric functions \& and $ satisfy the same second-order dif- 
ferential equation, the general solution of Eq. (|3.6j ) is 



^(i) = t [C $(1 + |; ±*t 2 ) + C" (1 + \ia 2 2l |; |zi 2 )] . (3.10) 
This is one convenient form for the solution for iV = 2. 

4 Case N = 2 — an Alternative Approach 

In the existing treatment in the literature for linear electron density and two types of neu- 
trinos [I], as reviewed in Sec. 3, the crucial step is to recognize that the second-order 
differential equation Q3.4|) can be solved exactly in terms of known higher transcendental 
functions, either parabolic cylinder functions or confluent hypergeometric functions. More 
generally, for N types of neutrinos, the corresponding differential equation is of iVth order. 
Even for iV = 3, the third-order differential equation is not one for any well-known tran- 
scendental function. Therefore, in order to be able to generalize the treatment of N = 2 to 
larger values of N, we must recast the solution of Sec. 3 so that parabolic cylinder functions 
and confluent hypergeometric functions do not play an essential role. 

A useful question to ask is the following: In what way is the linear electron density 
especially simple? The answer must be sought in Eq. Q2.6p , from which it is seen that the 
independent variable t appears only in one matrix element, and furthermore, it appears only 
linearly in that element. This implies that, if Fourier transform is applied to the differential 
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equation ( |2.4j ), the differentiation with respect to the Fourier-transform variable appears 
only once. Hence it is expected that an explicit expression can be obtained for the Fourier 
transform of ip. 
Let 

F(Q = 7T dte^Mt), (4.1) 



2tt _ 

then it follows from Eq. (|3.4j ) that F(() satisfies the first-order differential equation 

-C 2 ^(C) - |hcf(()] + (4 - i)F(0 = o, 

where we have omitted all terms from t = ±oo. This differential equation simplifies 
immediately to 

»'C - (C 2 - 4)F(0 = o, 



or 



1 dF(C) i, 2 , 2 



F(C) d( ( 



4-C). (4.2) 



Integration over ( gives 



F(C) = const. e -^/ 2 ( ia l (4.3) 



From the inequality ( |2.11| ) , it is seen that the function on the right-hand side of Eq. 
has a singularity at 

C = = b 2 . (4.4) 

Therefore the constant in (|4.3|) can take on different values for £ positive and for ( negative. 
In other words, the differential equation ( |4.2| ) is really two differential equations, one for 
£ > and the other for ( < 0, consistent with the fact that the right-hand side of Eq. ( |4.2| ) 
has a singularity at £ = 0. With this observation, it is natural to define 

-iC 2 /2 (_/-\ial f or /" < n 

F(0 = C { Q f* (4.5a) 

for C > 0, 



and 
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Inverting the Fourier transform ( |4.1| ), this choice leads to 

r-0 



c/Ce" iC *e-^ 2/2 (-C) iai 



— oo 
oo 



(4.6) 



With the notation ( [2.12| ), these two formulas ( |4.6| ) can be written as 



-Kt e -< 2 /2 |£| ia l ? 



(4.7) 



for n — 1,2. 

It remains to show that both x/j^ (t) and ip^ (t) are confluent hypergeometric functions 
of the correct parameters and argument. For this purpose, it is convenient to define 



so that it follows from Eqs. flOf ) that 

^\t) = tfj c (t)+t^ s (t), 

^f\t) = ^ c (t)-#.(t). 



/■oo 

/ dCco S ((t)e-* 2 / 2 e a * 
Jo 

d( sin(Ct)e- ic2/2 C ai 



It is found that 



V> c (t) = e— / 4 e™^ 4 2(~ 1+ ^)/2 r( | + | ia 2) $( I + i m 2 ; i. i^ 2) 



and 



^(t) = _ie^/ 4 2^/ 2 r(l + M)t$(l + i 



,2 3. 1 - + 2 



2 m 2i 2' 2 



2t 2 ). 



(4.9) 



(4.10) 



(4.11) 



There are various ways to verify Eqs. (|4 . 1 0|) and (|4.11| ), including carrying out power series 
expansions in t for the left-hand and right-hand sides. 

Finally, we note from Eq. (7) on p. 257 of reference || that 

(a, c; x) = ^^~ C \. $(a, c; x) + ^ x x ~ c $(a - c + 1, 2 - c; x). (4.12) 



r(a-c+l) 



r a 



Therefore, the results of Sec. 3 and this section are the same. 
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5 General values of N 



The procedure of Sec. 4 for N = 2 can be generalized in a straightforward way to larger 
values of N. Indeed, this is the major advantage over the previously known ones as reviewed 
in Sec. 3. This generalization to arbitrary values of N is to be carried out in this section. 
Thus, the differential equations ( |2.4|) need to be solved under the constraints (p. 71 ), (|2.10| ), 
and (grip . 



By Eqs. (|2.5| ) and (|2.6|) , the Eqs. (|2.4! ) are more explicitly 



#i(*) 
dt 



N 



i=2 



and, for k = 2,3,4... N, 



(5.1) 



In order to get a differential equation for ipi(t), apply the operator 

N 



n 

k=2 



■ d h 
l dt~ hk 



to Eq. flOl) . By Eq. this gives 

N 

n 



k=2 



■ d h 
I— Oi. 

dt h 



d 



x N N 

/ 4—0 L—O 



dt 



j=2 k=2 



i ~dt ~ bk ' 



(5.2) 



(5.3) 



Eq. ( |5.3| ) is an iVth-order ordinary differential equation for ipi(t); it reduces to Eq. 
when N = 2. 

Following Sec. 4, define the Fourier transform F(() of ipi{t) by Eq. (|4.1|) , then the 
first-order differential equation for F(() is 



N 



k=2 



f ft \ N N 

c-^ U(o = E a iIl(c-wc) 

\ ' ' 4 — I O 



(5.4) 



j=2 k=2 
k+0 



or 



(A \ N 2 



(5.5) 



or 



1 dF(C) 

no <k 



v «l 



i=2 



(5.6) 



S 



This Eq. ( |5.6| ) is the generalization of the previous Eq. ( |4.2| ) for N = 2. Integration over ( 
gives the generalization of Eq. ( [4.3|) : 



N 



F(0 = const, e 



X 2 /2 



II(C - K 

3=2 



(5.7) 



From the function on the right-hand side of this Eq. ( |5.7|) has singularities at 



for j — 2,3, ... N. Therefore, define for n — 1,2,3 ... N 



F { l ) = { e"< 2 /2 n ^ =2 | C _ hr ] for bn<(< hn+l) 
otherwise. 



(5.9) 



Inverting the Fourier transform fli.lD then gives the desired N linearly independent solu- 
tions of the differential equation ( |5.3| ) as 



j-bn+X 



N 



i=2 



(5.10) 



for n = 1,2,3, ... N. In both Eq. O and Eq. ( pT0|) , the notation of ( g3g ) has been 
used. The general solution of ( |5.3|) is of course 



JV 



^(t) = X) cr »^ B) (t) 



(5.11) 



n=l 



with arbitrary constants C n . 

The other components of the ^(t) of ( |2.5|) can be easily obtained also, and the result is 



where 



N 



^(t) = ^C^ (n) W, 



n=l 



(5.12) 



rbn+i 

i) {n \t) = I d(e~ Kt e 



-Kt P -K 2 li 



N 

n K - 

\j=2 



1 

C~b2 
C-&3 



a N-l 

L (-b N 



(5.13) 
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6 Limiting Behaviors for Large Distances 



The next task is to obtain the limiting behaviors of the various components of the wave 
function when the distance x is large, either positive or negative. In other words, the 
problem is to find the limiting behaviors of the ipj n \t) } as given explicitly by Eq. ( |5.13| ), 
both for t — > — oo and t — > oo, with all the a's and b's fixed. It is important to remember 
that these two limits correspond to different physical problems, as discussed after Eq. (|2.3|) . 

The consideration here will be limited to the part of the asymptotic behavior that does 
not vanish as t — > ±00. This is the physically interesting part. There are two possible 
types of contributions, from points of stationary phase and from end points of integration. 

6.1 Points of Stationary Phase 

From Eq. ( |5.13| ), the points of stationary phase are determined by 

|(-a-ic 2 ) = o (6.1) 

or 

C = -t. (6.2) 
In Eq. (|6.1J ), the additional phase due to the factor 



N 



UK 

i=2 



-bj\ 



is not included because the aj and bj are all fixed while t —>■ ±00. Eq. ( |6.2| ) implies that 
this point of stationary point is relevant only to: 

• ip^it) as t — > 00, and 

• ip[ N \t) as t — > —00 

in view of Eq. ( |5.13D . For example, when j > 1, ipj(t) as t —>■ 00 and ipj N \t) as t —>■ —00 
both behave as 1/t in absolute value so far as the contribution from this point of stationary 



phase (6.2) is concerned. 



6.2 End Points of Integration 

It is seen from Eq. ( |5.13| ) that, when k > 2, there is an extra factor of 



associated with ^ (t) . But the range of integration for this ip^ (t) as given by Eq. ( |5.13 ) 
is from b n to b n+ \. Therefore, the contribution from these end points of integration can 
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lead to a non-zero answer only when the index k appearing in the expression ( |6.3| ) agrees 
with either norn + 1. In other words, there are non-zero contributions as t — > ±00 only to 
ip k k ^(t) [i.e., n — k — 1] and [i.e., n = These particular components are given 

by 



■k-1 



N 



bk+l 

h 



II(c - b j) ia2j U 1 ' - o ia 

L i=2 - 1 L j=fc 

r fc n r W 

I1(C - hY a > II - 4 

i=2 - 1 L i=fc+i 



-Qfc 



(6.4) 



These Eq. ( p.4| ) are exact. 

Since the important contributions come from the vicinity of £ = all the £'s in 
Eq. ( |6.4j ) can be replaced approximately by b k except in the factors e _i ^', bk — C, and ( — bk- 
Therefore 



2V 



a 2 
Jk\ 3 



-ib 2 k /2 



■j=2 

N 



(-a k ) I dCe-^ih-Q 



-l+iat 



j=2 



(6.5) 



or 



N 



a 2 

>k\ 3 



L j=2 

N 



-ib 2 k /2 -ib k t 



n\ bj -b 



2 



■j=2 
3& 



-a k ) I dxe iXt x~ l+l < 



ixt — 1+iar 



a k I dxe x 




(6.6) 



This integral can be evaluated exactly in terms of the gamma function. 



6.3 Results 



Fig. 1 shows which ones of the various ip^ft) have non- vanishing behaviors for t 
and t — > 00. 

These non-vanishing behaviors are: 
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Figure 1: Table of non- vanishing components of ^ (t) as t — > ±oo. A cross means that 
the component is non-vanishing both for t — > — oo and t — > +oo; the symbol oo means for 
t — > +oo only, and — oo means for t — ► — oo only 



For £ positive and large, 



A' 



i=2 

■ A 



■3=2 
3+k 



while, for t negative and large, 



A' 



V£°(f) 



-ibl/2 -ib k t 
~ 6 fc 6 



•i=2 

• A 



•i=2 



n - 6 * 



(-a fe ) e -^/ 2 r(i^)r^, 



n i & i - ^ 



n \ b i - ^ 



(-a fc )e^/ 2 r(za 2 fc )|t|- ia S 



a k e-™l/ 2 r(ia 2 k )\t\- ia l. 



(6.7) 
(6.8) 

(6.9) 



(6.10) 
(6.11) 

(6.12) 



All the other components approach zero as t —>■ oo and as t — > — oo. In the asymptotic 
formulas ( |6T7| ) and ( |6.1U[ ), a is the quantity 

A 

a = ^a). (6.13) 



i=2 
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In the formulas (|OD , (IBTTT 



) and ( |Q2|) , 



2 < k < N, (6.14) 
where iV as always is the number of neutrino species. 

7 Discussion 

When we started to investigate the MSW differential equations for three neutrino species 
in the case of the linear electron density, we were mostly interested in various possibilities 
of finding approximate solutions. Therefore, it was quite a surprise to us that these coupled 
differential equations can be solved exactly not only for three, but also for any number of 
neutrino species. 

In the work of Wolfenstein, Mikheyev, Smirnov [|l] and others on the sun taking into 
account two species of neutrinos, it has been found that most of the effect takes place in 
a fairly narrow region around a particular value of the electron density. Because of this, it 
is quite accurate to use a linear approximation to the electron density. 

For more than two species of neutrinos, it is no longer true in general that there is a 
narrow region for most of the activity. Nevertheless, there are a number of circumstances 
where this is true. However, the conditions for this to hold has not yet been studied 
systematically. This is one direction for future work. 

Under the assumption of the electron density being a linear function of distance, the 
exact, general solution of the MSW differential equation is given by Eqs. (|5.12|) and (|5.13| ). 



This solution is in the form of a number of single integrals. When the number of neutrino 
species is more than 2, these integrals cannot be evaluated in terms of known functions, 
and therefore their properties need to be investigated. A small step in this direction has 
been taken in Sec. 6, where the asymptotic behaviors of these integrals have been evaluated 
for large distances but with all the a's and o's held fixed. It is believed that, in so far as 
this case of linear electron density is applicable to the physically interesting case of solar 
neutrinos, the asymptotic evaluation of Sec. 6 is far from being sufficient. It is more likely 
that not only the distance, but also some of the parameters, the a's and 6's, are large. This 
is a second direction for future work. 
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